Inter-species tunneling in One-dimensional Bose mixtures 
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We study the ground-state properties and quantum dynamics of few-boson mixtures with strong inter-species 
repulsion in one-dimensional traps. If one species localizes at the center, e.g., due to a very large mass compared 
to the other component, it represents an effective barrier for the latter and the system can be mapped onto 
identical bosons in a double well. For weaker localization, the barrier atoms begin to respond to the light 
component, leading to an induced attraction between the mobile atoms that may even outweigh their bare intra- 
species repulsion. To explain the resulting effects, we derive an effective Hubbard model for the lighter species 
accounting for the backaction of the barrier in correction terms to the lattice parameters. Also the tunneling 
is drastically affected: Varying the degree of localization of the "barrier" atoms, the dynamics of intrinsically 
noninteracting bosons can change from Rabi oscillations to effective pair tunneling. For identical fermions (or 
fermionized bosons) this leads to the tunneling of attractively bound pairs. 

PACS numbers: 03.75.Mn, 03.65.Xp, 03.75.Lm 



I. INTRODUCTION 



Ultracold atoms have become an important tool in the in- 
vestigation of strongly correlated quantum systems |3]. Be- 
sides controlling the external degrees of freedom via opti- 
cal and magnetic traps, one can vary the atoms' interaction 
strengths by the use of Feshbach resonances 0]. In one di- 
mension, confinement-induced resonances enable one to 
uncover the whole range of interactions from strong attrac- 
tion to hard-core repulsion. With this flexible toolbox at hand, 
it has become possible to study both the steady state and the 
quantum dynamics of atoms at an unpreceded level of accu- 
racy. For example, optical lattices have proven to be power- 
ful quantum simulators, giving insight into well-known con- 
densed matter phenomena, like, e.g., the superfluid to Mott- 
insulator transition iSHt]. While many studies concentrate on 
the weakly-interacting regime, in particular, in one dimen- 
sion, the strongly-interacting limit yields novel effects and 
valuable insights into properties of strongly-correlated sys- 
tems. If single-componentbosons in one dimension repel each 
other infinitely strongly, the bosons' hard-core repulsion can 
be mimicked by Pauli's exclusion principle and accordingly, 
they can be mapped to noninteracting fermions jg]. While the 
hard-core bosons are identical with the fermions in all local 
quantities, like the density or the energy spectrum, nonlocal 
properties such as the momentum distribution are very differ- 
ent. This phenomenon, known as the Tonks-Girardeau gas, 
has attracted a lot of attention [lfl, especially after its experi- 
mental realization 

More recently, one-dimensional bosonic mixtures - consisting 
of different species or hyperfine components - have received 
much attention [M [IlL E2, SI, [M Ell [TE E3, El] (see also 
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Refs. therein.) Their rich physical properties are rooted in the 
interplay between intra- and inter-species forces. Moreover, 
the trapping potential, in particular the frequencies of the con- 
fining potentials, may be made species dependent. While the 
research focus has been on static properties in most investiga- 
tions, it would be thrilling to find out more about their quan- 
tum dynamics. More specifically, we want to study a novel 
kind of tunneling lfl9ll in this paper. Most setups assume that 
the atoms tunnel throu gh a classical barrier provided by an 
external potential (see El |H HI HI HI IH] and Ref. 
therein). Only recently, studies have put the spotlight on dy- 
namics through quantum systems, such as atomic quantum 
dots ll27l l28ll and quantum barriers ll29l l30ll in optical cavi- 
ties. Here, we investigate tunneling through a material barrier 
realized via a quasi-one-dimensional mixture of two atomic 
species, where each component experiences a different degree 
of localization. If one species is tightly confined to the cen- 
ter of the trap, we show that for strong inter-species repulsion 
it can indeed be understood as an effective tunnel barrier for 
the other one, revealing the same properties as a system with 
a classical barrier. As that confinement is relaxed, the "bar- 
rier" atoms move due to the backaction of the other species. 
The dramatic effect of this additional correlation on the mo- 
bile component is studied both for the ground state and for the 
inter-species tunneling dynamics. 



Our paper is organized as follows: Sec. [TT] introduces the 
setup. The ground-state properties of the system are discussed 
in Sec. [Till first for the static case, where one species is com- 
pletely frozen out at the center of the trap (Sec. MI St . then al- 
lowing for some flexibility of the "barrier" atoms (Sec. lIIICl ), 
completed by an explanation of the observed effects within 
effective lattice models. Section HVl focuses on the quantum 
dynamics of two identical atoms interacting with one "bar- 
rier" atom, rounded off by the discussion of higher particle 
numbers in Sec. lIV Dl 
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n. SETUP 

We consider a mixture of two distinguishable bosonic 
species, which we will label a = A,B. These may cor- 
respond to different hyperfine states of one and the same 
species in the case of identical masses, or to two differ- 
ent atomic species or isotopes in which case the two con- 
stituents could also have different masses. We study quasi- 
one-dimensional systems, where we assume the atoms to be 
tightly confined in the transverse direction, such that these de- 
grees of freedom can be integrated out. The Hamiltonian for 
an arbitrary mixture of N = Na + Nb bosons then reads 
H = J2a=A b H a + ^ab with the effective one-dimensional 
single-species and inter-species Hamiltonians 
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The effective one-dimensional intra- and inter-species interac- 
tions (see JU for more details) can be written as 
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where the transverse trapping frequency ui± is assumed 
to be the same for both species. The 3D scattering 
lengths are denoted as Oo.ab and ao.o-, while a±. a = 
y // 2h/m a ujj_; a± t AB = \fhj Hab^±_ describes the transverse 
confinement length, pab — toawb/(toa + ^b) is the re- 
duced mass and C « 1.4603 a constant. 

In the following we will focus on harmonic trapping po- 
tentials U a = ^m a uj^x'^ and repulsive forces per, 5ab > 0. 
Rescaling to harmonic oscillator units a a = yJhjmA^A of 
the "A" component enables us to eliminate tua = ua = h = 
1 by exploiting the scaling 



Hoirna^cg^Xo) = frhJAHo ( — ^, — ,g'„\X' a 



with X' a = X a /a A and g' a = g a ^m A /h 3 ujA, g' AB = 
gAB \J wa/^wa. This leaves us with a Hamiltonian free of 
any dimensionful parameters, which is very convenient for nu- 
merical simulations. In terms of the mass ratio a — ,J1A - and 
the frequency ratio (3 = the rescaled Hamiltonian now 
reads 



i=l 

-9A Yl 

i<j;a=A, 



2af3 2 
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where all primes have been left out for convenience. For sim- 
plicity, let us assume that both species experience the same 



trapping potential {[3 = 1) while their masses may vary. This 
restriction is not crucial, and experimentally large differences 
in the trapping potential (J3 « 0) might be more feasible than 
large mass differences (a « 0). 



III. GROUND-STATE PROPERTIES 

In order to investigate this two-component Bose mixture, 
we load both species into a harmonic trap of a frequency u>. If 
the masses of both constituents are equal, ijia = tub, the dis- 
tribution of the two species throughout the trap depends on the 
atom number of each constituent N a and the intra- and inter- 
species interactions g a , gAB- In addition, a variable degree of 
localization for each species can be realized by changing the 
mass ratio a. A large mass difference toa -C tub (as, e.g., for 
Li-Cs) leads to different confinement lengths for both species 
such that <2b 1. This results in a separation of length scales 
between the lighter A and the heavier B bosons, such that 
each sort only experiences an effective potential aroused by 
the other one. In particular, the case of strong inter-species re- 
pulsion, gAB 3> 1, is interesting, and we will concentrate on 
this scenario throughout this paper. In that case, the B bosons 
act as an effective barrier for the lighter A component. While 
this barrier remains localized at the center of the trap for very 
unequal masses (a f» 0), the B bosons begin to move due to 
the backaction of the lighter species for larger ratios a > 0. 
We will show how the density profiles change for a varying 
confinement length starting from equal masses, a = 1, in 
Sec. IIII Al Thereafter, we will describe the limit of a static 
barrier, a — > 0, in Sec. IIII Bl and carry on with an investiga- 
tion of the effects of varying mass ratios a > in Sec. IIII Cl 



A. Confinement-induced Demixing 

Let us begin with the completely symmetric case where 
both species are identical in masses, external trapping and 
atom numbers (a = j3 = 1,Na = Nb), which may corre- 
spond to two different internal states of the same species (e.g. 
two different hyperfine states). This setup is depicted for 
Na = Nb = 2 bosons in Fig. Q] Accordingly, the Hamil- 
tonian H reveals complete symmetry between the bosons A 
and B, resulting in identical density profiles: for both species 
these are spread out over the entire trap with a dip at the ori- 
gin x = resulting from the strong inter-species repulsion 
gAB = 25. This statement has to be taken with a bit of cau- 
tion: the densities pa(x) = Pb{x) are already ensemble av- 
erages: due to the strong inter-species repulsion, in a single 
measurement, we will always find all atoms Na at one side 
of the trap and all atoms TVb at the other side, which is de- 
picted in the two-particle density p A ^(x\, x%) in Fig.Q] The 
probability of finding an A and a B boson simultaneously at 
the same side of the trap consequently converges to zero as 
gAB — * oo, and one may think of this state as a superposition 
of the form * = \N A , 0) ® |0, 7V B ) + |0, N A ) ® \N B , 0) 0, 
where |JVr, Nl) denote a state with Nb, atoms at the right and 
Nb at the left side of the trap. 
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Figure 1: (color online; dimensionless harmonic oscillator units a a 
throughout) Top: Influence of the mass ratio a on a system with 
g a = 0, qab = 25, Na = Nb = 2. One-particle densities of the A 
bosons (left) and the B bosons (right) for a varying mass ratio a (be 
aware of the different scales in the plots). Bottom: Two-body density 
Pab( Xi > X2 ) f° r Q = L 0.5, 0.1 (from left to right) 



This symmetry is easily broken and the hidden phase sepa- 
ration becomes obvious as soon as we change the relative con- 
finement by varying the mass ratio a. Let us focus on the latter 
as depicted in Fig. Q~| for any q ^ 1 the symmetry is broken 
and the heavier B bosons experience a stronger confinement. 
Their oscillator length is reduced by a factor of as = \[u, 
preventing them from moving to the sides of the trap. For a 
decreasing mass ratio a, they consequently begin to assemble 
more and more at the center of the trap, sandwiched between 
the lighter A bosons that move to the sides. Accordingly, the 
dip in the density of the B bosons becomes less pronounced 
as a is decreased: while it is still visible for a = 0.5, it has 
vanished for a = 0.1, and the density profile of the B bosons 
begins to resemble a narrow Gaussian, making clear why it is 
appropriate to treat them as an effective barrier. On the other 
hand, the dip in the distribution of A bosons becomes more 
pronounced for a decreasing a as a consequence of their ten- 
dency to avoid the B bosons. 



B. The static limit a w 

Imagine now the situation where one species (say B) is 
much more strongly localized at the center of the trap. As 
mentioned before, let us assume much heavier B atoms but 
equal frequencies, resulting in very different length scales 
for the two species. Consequently, the B bosons are effec- 
tively frozen at the center, while the density of the A bosons 
changes slowly throughout the entire trap. To get a better un- 
derstanding of this situation, let us now account for the B 
bosons within an effective model. Due to the different masses 
(a w 0), it is possible to carry out a Born-Oppenheimer ap- 
proximation and use a product ansatz for the wave functions, 

such that the densities can be written as jr^ f=a P~a A ^ ® Pb^ ■ 
We are now able to integrate out the B component to obtain 
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Figure 2: (color online) In the limit of a vanishing mass ratio a — 
0.001, the lighter Na = 4 bosons develop from an uncorrected 

profile ( — ), over a widened distribution ( ) to a fermionized state 

(■ ■ ■ ) for an increasing intra-species interaction qa and strong inter- 
species repulsion qab = 25 , while the density of the B bosons 
represents an effective barrier for the lighter species (- - -) that barely 
changes for different qa- 



an effective Hamiltonian for the lighter species 

H ( ° ] =H A + ti-B [H A bP { b b) ] = H a + sabX/^a,*)- (1) 

i 

This effective single-particle potential for the A bosons con- 
sists of the harmonic trapping potential and a sharp barrier at 
x = of width as given by the one-body density pb(x) = 
Ub(x) /Nb- Its height is proportional to the inter-species cou- 
pling and for strong repulsion <?ab ^S> 1, the system bears 
strong resemblance to an external double well potential. We 
will now proceed and compare this analytical model to our 
numerical results in the following subsections. 

1. Density profiles 

Let us first consider the one-body densities p(x), which de- 
scribe the probability of finding one particle at position x. The 
mass ratio has been chosen as a = 0.001, giving the density 
of the heavier B boson a <5-like shape: they simply experience 
a harmonic oscillator potential, with the ground state being a 
normalized Gaussian of width (with h = 1 and u> = 1): 



Pb(xb) = 




converging towards 5(xb) for a — » 0. The inter-species inter- 
action strength has been fixed at ^ab = 25, while the intra- 
species interaction #a is varied. This barely affects the heavy 
B bosons, which remain unchanged to a good approximation 
on variations of gx- 

On the other hand, the strong inter-species repulsion ren- 
ders the central region energetically forbidden for the lighter 
A bosons. From their perspective, the harmonic trap is con- 
sequently split into two halves and well described by the ef- 
fective model (eq. (Q~|i), as we will witness in the following. In 
this sense, for g A — all bosons reside in the single-particle 
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Figure 3: Two-body density p AA (xi, X2) of the lighter A species 
for: Na = 4, Nb = 1 bosons, a mass ratio a — 0.001 with fixed 
Pab = 25 and increasing g\ = 0; 0.5; 25 (from left to right). 



ground state 'I' a = 4>t A of the effective model, and the den- 
sity distribution is consequently delocalized over the entire 
trap as illustrated in Fig. [2] If we increase qa to cja ~ 2, the 
distribution is broadened due to the intra-species repulsion. 
For even higher interaction strengths, a completely different 
picture is observed: a structure of Na/2 peaks on both sides 
of the effective double well emerges. Physically, this means 
that if we measure the position of an A boson, it is very likely 
to find it at Na discrete spots - and less likely to detect it 
anywhere in between. This localization has a plausible physi- 
cal explanation - due to their strong repulsion, the bosons try 
to isolate from each other. In free space, they would simply 
start to move apart, until they don't feel their repulsion any- 
more. In our setup they cannot do so, as they are confined 
in a harmonic trap with an effective central barrier consist- 
ing of bosons of a different species. Therefore they try to 
minimize their mean interaction energy and are consequently 
pinpointed at more or less discrete positions. In compari- 
son, the same phenomenon would arise for Na noninteracting 
fermions: due to the Pauli principle they are not allowed to 
occupy the same state and begin to localize in a configuration 
with a maximized mean distance. This behavior relies on the 
one-to-one mapping between hard-core bosons and noninter- 
acting fermions in one-dimensional systems that has been rig- 
orously proven in Girardeau's seminal paper in 1960 Jo and 
is commonly referred to as fermionization (see, e.g., Bill for 
more details). It relies on the fact that bosons with infinitely 
strong repulsion become impenetrable, which means mathe- 
matically that configuration space becomes disconnected into 
regions {xi 7^ Xj\i < j}, allowing us to apply the Bose-Fermi 
map ty + = with A = Yii<j s S n ( x i — x j)- K is both 

valid for stationary and explicitly time-dependent states and 
shows that hard-core bosons and noninteracting fermions will 
agree in their spatial densities and energies, while nonlocal 
properties like the momentum distribution are different. 



2. Two-body correlations 

To get better insight into the mechanism underlying the 
fermionization crossover, we also examine the two-body den- 
sities. In the infinite-mass limit a — ► 0, pa and p-Q decou- 
ple, and we concentrate on the two-body correlation function 

(2) 

Paa( x 1i x V- ft describes the probability of finding one par- 
ticle at position x\ and the other one at x-x. For qa = 0, 
correlations are absent and p AA — p A ® p^ . This is il- 



lustrated in the left plot in Fig. [3] where the depleted cross- 
like centered region results from the effective potential caused 
by the B species. Besides, the A bosons are coherently dis- 
tributed over both "wells". This changes as we increase qa to 
cja = 0.5: a correlation hole at the diagonal xa.i = xa,2 be- 
gins to form. It has the intuitive explanation that two bosons 
avoid being in the same place, as this would be very costly 
in energy. This phenomenon is only visible in the two-body 
density, while it is smoothed out in the one-body density 

Pa\ x ) = / dx'pAA{ x -> x ') by the integration over the entire 
space. 

Tuning the interaction strength up to the intra-species 
fermionization limit <ja = 25, we witness the same picture as 
in the one-body density: the bosons are localized at discrete 
spots and the two-body density forms a checkerboard pattern. 
The interpretation is as follows: if one boson is fixed at a dis- 
crete spot xa,i, the probability of finding a second boson at 
X A,2 = %A,i is zero, while finding it at three other discrete 
positions is equally likely. The same structure emerges for 
noninteracting fermions, again manifesting the equivalence 
between the boundary condition resulting from the strong re- 
pulsion of hard core bosons, *f?\ XA i=XA j = 0, and the Pauli 
principle for fermions. 



C. Beyond the static limit 

So far, we have assumed that the B bosons are heavy 
enough to be completely frozen out at the center of the trap, 
and it has been shown that for a = ttia/tib ~ this 
"quantum barrier" leads to the same ground-state properties 
as an external double well. For less restrictive mass ratios, 
the B atoms are expected to move due to the backaction of 
the lighter A atoms. They try to avoid the strongly repul- 
sive lighter atoms and move towards the regions in the trap 
where pa is smaller. Thinking of the system as an effective 
double well, this means that the well, in which the A bosons 
reside, is now wider than the empty well. Hence, this ren- 
ders configurations where all A atoms reside at one side of 
the trap energetically more favorable than a setup, where the 
A bosons are distributed equally over both sides. This addi- 
tional correlation will have an immense effect on the ground 
state properties of the A bosons. We want to investigate it in 
more detail and we will therefore first derive effective models 
valid for both weak and strong intra-species interactions c/a in 
Sec. lIIICTI and compare it to our numerical findings thereafter 
in Sec. |IIIC2| 



1. The effective model 

The effective Hamiltonian 
The expansion of the Hamiltonian in orders of a is highly 
nontrivial as H(a) by itself does not suggest a straightfor- 
ward power series in a. In order to write the Hamiltonian in 
a form that is suitable for this expansion, we introduce center- 
of-mass and Jacobian coordinates relative to Xb- Without loss 
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of generality, we set N B = 1 and transform to the new coor- 
dinate system (r± , . . . , r NA _ i , R) 



+ 0(a 3 ), 
V A (x 1 ,x 2 ) 



which 



adds 

- x 2 ). 



to the local interaction 



Yh=1 m A,i x A,i + m B x B 



n = XA,i - x B ; r k = XA,k 



i=l m A,- 



m B 



SiI A i m AXA,i + m B x B 



i?= 



with the corresponding reduced and total masses 
Yh-i m A,i + m B ) m k 

m& = — =^ ; M = Y, l= i to a, ; 

2^=1 m A,i + m B 
The relative Hamiltonian thus reads 



m B . 



(A. 



i 



5AB $^(cfc,i fa}). (3) 
fc<i 



where the term containing the intra-species interaction has 
been left out for simplicity for now. The last term in eq. (01 
consists of combinations of relative coordinates of the form: 



SAB 2^ 

k<l 



S(c k i {r k }) « .9ab(<5(?"i) + <5(r 2 + an) 

+S(r 3 +an + ar 2 ) + ...), (4) 



containing the perturbation parameter a. A Taylor expansion 
of eq. (01 leads to a Hamiltonian in orders of a 



with H<® = E fe (£ + WD + 5ab E fe *fa). 



rl\ = Efe<i3A B ar fc( 5'(r ; ), 



(2) 
rcl 



J2k<i ^9ABOt 2 rk5"(ri), etc., where the prime de- 
notes a derivative dS(rk)/drk- Assuming that rj, « x\,k~x B 
for small a, we can trace out the B bosons to obtain an effec- 
tive Hamiltonian for the A bosons: 



H A (a) 



H 



(o) 



O: 



2 9AB i( x A,k ~ x A,i )n B (x A: i) 

k^l 



- n B {xA,i)} + ~~~~ $AB y^j2n' B (xA,i) 

k=jtl 

+ (x A ,i - XA,k)n B (x A ,i)} + 0(a 3 ). 



(5) 



The permutation symmetry between the A bosons has been 
broken by the assumption rk ~ x A ,k — x B and has con- 
sequently been restored by hand in the above Hamiltonian. 
While the first term in eq. (0 recovers the initial infinite-mass 
approximation, the higher order terms may be under- 
stood as an additional external potential on A: SUa(x) = 

Sab^V 1 [~a[xn' B (x)+n B (x)} + f [2n' B (x) + xn B {x)\ } 

+ 0(a 3 ) plus an induced nonlocal interac- 
tion between two A atoms: 8Va{x\,X2) 



An effective Bose-Hubbard model 

Given the multi-well (here: two sites) geometry experi- 
enced by the A atoms, the use of simplified lattice mod- 
els suggests itself. For sufficiently low intra-species interac- 
tions g A < 1, an effective Bose-Hubbard model proves to 
be useful. It relies on the lowest-band (limiting its use to 
small interaction strengths) and the tight-binding (limiting the 
model to deep wells) approximation, which translate in our 
setup into a relatively small intra-species interaction strength 
gA and a sufficiently large inter-species interaction strength 
<?ab- The eigenstates of the effective one-body Hamiltonian 
h A = — \d\-\-\x\-\- gABn B (xA), describing the completely 
localized case a — > 0, readily give access to a basis for this 
effective Hubbard model. These eigenstates are arranged in 
bands (in analogy to the Bloch functions for lattices) here de- 
noted by (3, consisting of doublets of symmetric (x) and 

antisymmetric (f>\ (x) states. These delocalized functions are 
unsuitable for the evaluation of local properties and we apply 
a unitary transformation to obtain localized Wannier functions 

wi 0) (x) = (4>i\x) ± 4>f\ x )) instead, where s = ±1 
denotes the "lattice" site. Let us now expand the effective 
Hamiltonian (0 Ha{q) in its second-quantized form 

H A (a) = J dx¥(x)[h A + SUA(x)}y(x) 

+ \ J dx 1 dx 2 ^(x 1 )^{x 2 )[gA8(x 1 -x 2 ) 
+ SV A (x 1 ,x 2 )}%(x 2 )$(x 1 ) 

in terms of the Wannier functions $>(x) = Y^ s ®s w <* ( x )' 
where a s describes the annihilation operator at site s. We only 
take into account the lowest band as gA is sufficiently small 
and finally obtain the effective Bose-Hubbard Hamiltonian 



H { r = -j (o) e # 

(s,s>) 



]Tn s (ri s -l). (6) 



The first term is the hopping matrix element between nearest 
neighbors j(°>, which is renormalized by SUa, yielding 



J(°) = (wi [)) \h A +SU A (x)\w^) 



-J, 



(0) 



<5J(°) (and 



correspondingly for the other bands (3). While the correction 
term 5J(°\ tends to be small, the on-site renormalization 
due to SV A , u = (wi 0) ® 2 \V A + 6V A \wi 0) ® 2 ) = u Q + Su, 
can have a huge effect on the system: it induces correlations 
between different A bosons and effectively leads to an 
additional intra-species coupling. To make this more clear, 
we carry out some analytical estimations by assuming a 
Gaussian density profile given by eq. (0 for the B bosons 
and harmonic oscillator functions for the Wannier states 
with a displacement x s and a width 7. In this approxi- 
mation, we find that the correction to the tunnel coupling 
<5J(°) = -V2 5AB a 2 cxp(~x2/ 7 2 )/(5 7 3 ) + 0{a 3 ) is 



suppressed by a factor a 2 and consequently small. Besides, 
24^- ja [xin' B (x 2 ) + x 2 n' B (xi)] — ^- [xin B (x 2 ) + x 2 n B (xi)] > a \\ interaction terms involving different sites are several 
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orders of magnitude smaller than the on-site interaction 
and will hence be neglected in the following. By con- 
trast, the correction to the on-site interaction can have a 
remarkable effect on the system if it is of the same order 
as the a = contribution. The most remarkable feature 
is that all Su are smaller than zero for any choice of 7 and 
x s as it is visible in the contribution to first order in a, 
8u = -^=^ag AB x 2 s cxp(~x 2 J(a + 7 2 )) + 0(a 2 ) . This 
demonstrates that the induced correlation can be interpreted 
as an effective intra-species attraction, increasing as a 
becomes larger. 

An effective Multi-Band-Hubbard model 
The lowest-band approximation naturally breaks down as we 
go to higher interaction strengths. This becomes graphical in 
the fermionization limit gA — » 00, where the hard-core bosons 
can be mapped onto noninteracting fermions. According to 
the Pauli principle these are prevented from occupying the 
same band and consequently, N A bands have to be included 
for calculations with N A hard-core bosons. Consequently, we 
can write down an effective Fermi-Hubbard model (FHM) for 
identical fermions: 



H 



(FH) 



E ^/W^-Em 



1 

2 



s,0 



-E E ^ (Q ^ Q ' ' ) /i Q)t /i ,3)t / s ( ' 9 ' ) /i a ' ) , 

s a ,p,a'P' 

where Na bands contribute as the fermions obey Pauli's 
exclusion principle and p!f^ denotes the on-site ener- 
gies. We have assumed that only nearest-neighbor tun- 
neling gives a contribution and that inter-well interactions 

Su^ a 13 ' are negligible. The operators /j' 3 ^ ^/i^J cre- 
ate (annihilate) a particle in the /3th band at the sth 
site, and its tunnel coupling is given by -J0») = 
(wi^\h,A + SUa(x)\uu^}. The on-site interaction is given 
hySu (a(3a'p') = ( w (<x) w (P)\SVA\w( a '^w^">) and the total in- 
duced interaction between different bands is quite elaborate 
containing all different combinations of creation and annihi- 
lation operators of the TVa bands. In order to make this term 
more comprehensible we write it in a different form 



H int = \ E SvW^nW 

s,a</3 

+ i E SvW^W^fiWft 

Sv ( a 0a>'P') jja)f j(l3') j{ a ')^ 



E 



with the induced inter-band interaction strengths 

5v (a f 3a ' fS ' ) = (w (a) w^ ) \5V A \w {a ' ) w^' ) ) 
+ (w {m w {a) \SV A \w {(3 " ) w ia ' ) ) - (w^w^ISVaIw^w^) 

-{w^w^ISVaIw^w^}. 
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Figure 4: (color online) Beyond the static "barrier": Two-body den- 
sity p A l(xi_, X2) for N A = 4 and Nb = 1 bosons for (top) g A = 0, 
with increasing mass ratios a = 0.001, 0.002, 0.02 (from left to 
right); (bottom) g A = 0.5, a = 0.001, 0.12, 0.2. 



In this form, it becomes clear that in the case of two A 
bosons the interaction term can be reduced to i? on _ s i te = 

J2 S Sv^ owl ^hsha. In analogy to the bosonic case, we 
carry out some analytical estimations in a Wannier basis con- 
structed of harmonic oscillator functions. While the correc- 
tions to are small (see previous paragraph), the induced 
§ v ( a P a ) can become remarkable for sufficiently large a. 
The induced interactions are again negative, which means that 
the fermionized bosons effectively attract each other. The first 
order contribution between the lowest bands is of the form 
^ (0101) = g A Bax 2 s cM~^s/l 2 )(-2x 2 s + 7 2 )/(V^7 5 ) + 
0(a 2 ), which is negative for an appropriate choice of the lat- 
tice parameters, i.e. x s > j/y/2. Of course, this induced 
attraction does not violate the Pauli principle as the fermions 
still reside in different energy bands. 



2. Results 

The effective Hubbard models suggest corrections to the 
ground state resulting from the induced intra-species inter- 
actions Su. Let us now examine how these induced ef- 
fects manifest in the one- and two-body densities. While 
the inter-species interaction is fixed at a large repulsive value 
(Jab = 25, the intra-species interaction is varied. Let us 
first consider the limit with vanishing intra-species interac- 
tions 9a = 0. In particular, the two-body correlation function 
(2) 

p AA (xi, X2) gives some interesting insight into these systems 
(see Fig. |4|. For a as 0, there are no correlations between 
the A bosons and the two-body density looks exactly like in 
the case of an external barrier. This picture changes dramat- 
ically as we vary the a m requirement: already for rela- 
tively small mass ratios, correlations are induced, making it 
more favorable for the bosons to assemble at the same site 
and the population of the off-diagonal decreases. Already 
for values as small as a = 0.02, the A bosons cluster and 
the probability of finding them at different sites converges to 
zero. The effective Bose-Hubbard model in eq. (O is at hand 
with an explanation for this phenomenon: the induced attrac- 
tive interaction Su between the A bosons makes it energeti- 
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Figure 5: Evolution of a system with qa — 25, gAB = 25, Na — 
4, Nb = 1; 7bp: One-particle densities for increasing a; Le/r: of the 
A bosons, Right: of the B bosons (be aware of the different scales 
in the plots); Bottom: Two-body correlation function p A Y A (x\,X2) 
a = 0.9; 1; 1.1 (from left to right). 



cally more favorable to assemble several bosons at the same 
site. The resulting state can be thought of as a maximally 
entangled or "cat state" where the measurement of one bo- 
son determines the outcome of all preceding measurements: 
\-qj A ) = -±= (I-ZVajO) + |0,7Va)), as the entire ensemble as- 
sembles either on the left or the right side. 

For larger interaction strengths g A = 0.5, the system is in 
a "Mott-insulating"-type state for tiny mass ratio a = 0.001. 
Having witnessed the effect of induced attractions for <ja = 0, 
the question naturally arises whether it might be possible to 
compensate the repulsion and bring the system back into an 
uncorrected state. Considering the two-body correlations in 
Fig.[4](bottom), this picture is indeed confirmed: starting from 
an anticorrelated "insulating" ground state for a ss 0.001, the 
population of the diagonal gradually increases until we end up 
with a seemingly uncorrelated state for a « 0.2. This behav- 
ior is predicted by our effective Hubbard model: 8V A gives 
rise to an induced attraction Su « uq which can compensate 
or even outweigh the "real" intra-species interaction V A . 

Let us now move on to the strongly repulsive fermioniza- 
tion limit. This case is well characterized by the effective 
Fermi-Hubbard model in eq. (|7]i for small mass ratios a. It 
predicts induced correlations between fermions in different 
bands. In order to be measurable in the ground state, this in- 
duced interaction has to be of the same order as the band gap 
e (!) _ e (o)^ making it favorable for a second fermion to occupy 
the same well instead of the other empty well. The induced in- 
teraction is naturally tiny for small a as illustrated in the spec- 
trum in Fig. [7] and it only leads to corrections much smaller 
than the band gap. Therefore, no effect is visible in the ground 
state for relatively little mass ratios, and we will consequently 
look at the evolution of the states for larger a. One inter- 
esting question is: What happens to a mixture of two initially 
fermionized components (g A = 25) as their mass ratio a tends 
to one? For a = 1 with g a = g A B = 25, the two-component 
mixture is related to a single-component hard-core Bose gas, 
that in turn can be mapped to a single-component ideal Fermi 
gas via the Bose-Fermi mapping 

sunn 

[H [nt]. The 



density profile is therefore the same as for an ideal single- 
component Fermi gas (see Fig. [5}: for five hard-core bosons 
TV = N A + N-q = 5 humps emerge. This behavior also man- 
ifests itself in the two-body correlation function: if one bo- 
son of either species A or B is measured at a discrete spot, 
the probability of finding a second boson at the same place 
is zero and there are four equally likely discrete positions 
for the remaining bosons. This special case of equal masses 
and interactions is highly degenerate: there is a permutation 
degeneracy between A and B bosons, and moreover in the 
limit gAB — ► oo, the ground state wave function also degener- 
ates with the case, where both A and B are fermionic by the 
Bose-Fermi mapping. Due to the symmetry between the two 
species the correlation functions p A ^(xi, X2) = p AA {x\,X2) 
are identical lfl5ll . 

This degeneracy is lifted as we introduce a small difference 
between the A and B bosons, for example by varying the 
mass ratio a. Even for small deviations from a = 1, the 
symmetry in the Hamiltonian between A and B is broken. A 
slightly heavier B component (see a = 0.9 in Fig. [5]) already 
localizes at the center of the trap and the lighter A bosons 
move to the outer parts. This is due to the high intra-species 
repulsion between the A bosons, which overwhelms the 
potential energy: even if the density of the A bosons is higher 
(N A = 4) and it is thus more expensive for them to move to 
the sides of the trap than it is for the B boson, their strong 
intra-species repulsion forces them to do so, sandwiching 
the B boson at the center. Further decreasing the mass ratio 
a has no qualitative effect on the system: the B bosons 
remain localized at the center of the trap with a decreasing 
oscillator length ae = \/a, while the A bosons, preserving 
the checkerboard pattern, tend to move slightly to the inside 
of the trap (not shown). In the opposite case a > 1, the B 
boson, as the lighter species, moves to the sides of the trap, 
surrounding the A bosons, that are concentrated at the center 
forming three wiggles (see a = 1.1 in Fig. [5]). For the A 
atoms, the presence of the heavier B bosons at the sides of the 
trap amounts to hard-wall boundary conditions imposed via 
the strong repulsion between the two species. The harmonic 
potential has a minimum at x = 0, and the hard-core bosons 
of species A therefore localize around the central region and 
form N A — 1 = 3 wiggles. The fourth A boson is delocalized 
over the entire allowed region due to symmetry requirements. 



IV. INTER-SPECIES TUNNELING DYNAMICS 

Having investigated the equilibrium properties of a bosonic 
mixture where the heavier of two species acts as an effective 
barrier for the lighter one, it would now be interesting to learn 
about the nonequilibrium quantum dynamics of this setup. Do 
the lighter atoms tunnel through the heavier ones and how is 
this motion influenced by the backaction of the "barrier" for 
varying mass ratios? To tackle these questions, we prepare 
the barrier bosons B at the center of the trap while the lighter 
A atoms are entirely loaded into the left-hand side. This can 
be achieved by either displacing the trap center for the lighter 
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Figure 6: (color online) Inter-species tunneling dynamics of nonin- 
teracting (qa = 0) JVa = 2 bosons through a strongly repulsive B 
atom (<?ab = 8) : Top: Evolution of the population of the A bosons 
for increasing mass ratios a — 0.001, 0.01, 0.122 (from left to right); 
Middle: Same as above but for the B boson; Bottom: (a) Relative 
population of the right-hand side over time, PA,n(t)', (b) Probability 
f>A,2(t) of finding two A atoms on the same side for a = 0.001 ( — ), 
a = 0.01 (- - -), and a = 0.12 (• ■ • ) 



4.5 
4 
3.5 



7.5 
7 
6.5 
111 6 
I 5.5 
LU 5 
4.5 
4 
3.5 
3 



— 4M- 412 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



species or by blocking the right-hand side by shining blue- 
detuned laser light on it. At release time t = 0, the asym- 
metry in the potential of the lighter atoms is lifted and they 
begin to tunnel through species B, that react and in turn mod- 
ify the effective potential of the lighter bosons. To explore 
the tunneling dynamics of the A bosons through the barrier 
atoms, we fix the inter-species interaction strength at a large 
repulsive value <?ab = 8 and vary the intra-species interaction 
9 A > 0, covering the full fermionization crossover. 



A. The noninteracting limit qa = 

We first focus on the case of noninteracting A bosons 
c;a = 0, as illustrated in Fig. [6] For a vanishing mass ratio 
a = 0.001, the B bosons remain fixed at the center of the 
trap, while we observe simple Rabi-type oscillations between 
the right and the left side, illustrated by the evolution of the 
density of A bosons in time (top) and the percentage of atoms 
at the right side of the trap pa,r(0 == J °° Pk{x\ t)dx. The 
observed tunneling time scales become clearer if we consider 
the spectrum. In the limit of tiny mass ratios a = 0.001, it 
is appropriate to use the effective Hamiltonian in eq. (Q]) as a 
means to explain the occurring effects. It treats the B bosons 
as an effective barrier for the lighter A component. The result- 
ing effective potential for the A bosons thus resembles a dou- 
ble well, for which the lowest energy states for N\ bosons in 
the noninteracting limit, g A = 0, are given by distributing all 
atoms over the symmetric and antisymmetric single-particle 
orbitals of the lowest doublet (also see 12511 ). For Na bosons 



Figure 7: (color online) Spectrum for Na = 2 and Nb = 1 bosons 
with varying a; Top: gA = 0, Middle: gA = 0.5 Bottom: gA = 25; 
(?ab = 25 throughout 



this leads to N A + 1 energies E m = Eq + Eq + 2m 

with m = 0, ...,N A and 2J(°) = e { " ] - e ° ] being the en- 
ergy splitting of the first band. For vanishing interactions and 
TVa = 2 particles, only TVa + 1 = 3 states contribute and the 
population of the right-hand side can be easily computed as 

Pa,r(*) = (*(t)|6(*A)|*(t)> 

= i(l + c ci(* |e(a:A)|*i)cos(woii) 

+ cic 2 (*i|e(.T A )|*2)cos(cji 2 £)), (7) 

where c m = (\I r m |\I/(0)) and the (02) contribution vanishes 
as only opposite-parity states are coupled. As illustrated in 
Fig. |7J for <7a = the lowest levels are equidistant and con- 
sequently only a single-mode Rabi frequency c^oi = W12 — 
2 J (0) contributes. 

Another measure, that gives a better insight into the 
dynamical properties, is the pair probability pa, 2(0 = 
I{xi x 2 >o} AAA 0^1; £2; t) dxi dx2, which is the integral over 
the time-dependent two-body correlation function, giving the 
probability of finding two A bosons at the same side of the 
trap. Noninteracting bosons in a static double well potential 
tunnel independently as reflected in the two-body correlation 
function: if both atoms start out in one well, their pair proba- 
bility will drop from one to one half at equilibrium points of 
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the oscillation as illustrated in Fig. [6] 

Increasing the mass ratio, it becomes clear that even for 
small a = 0.01 the mobility of the B bosons has dramatic 
consequences for the tunneling dynamics of the lighter atoms: 
beats with a frequency W12 — woi begin to evolve on top of the 
sinusoidal oscillations. It is visible in the spectrum that a>oi 
is slightly smaller than U12 and eq. (|7]i explains the resulting 
two-mode oscillation. 

For even higher mass ratios a as 0.1, the tunneling slows 
down drastically with only a tiny faster modulation on top. 
The form of the spectrum (Fig. |7J gives an explanation for 
this behavior: for a = 0.12, the lower two lines Eo,i are 
virtually glued together to form a doublet, whereas the gap to 
E2 increases. This explains the occurring tunneling pattern: 
complete population transfer occurs at a fairly long time scale, 
determined by the splitting of the lower states woi, while W12 
only evokes a small oscillation on top of it. 

The effective Bose-Hubbard model (|6]l explains this phe- 
nomenon by introducing an effective attractive intra-species 
interaction Su < (see Sec. IIIIC lb . This attraction makes 
it more favorable for both bosons to assemble in the same 
well. The initial state, where all bosons are prepared in the 
left-hand well, is thus only in resonance with a configuration 
with all bosons in the right-hand well and not resonant with 
the state where one boson occupies each well. As the bosons 
cannot gain energy, they cannot tunnel to this energetically 
higher state, and single-atom tunneling is highly suppressed. 
In fact, the Bose-Hubbard model suggests the inte rpret ation 
of the bosons as attractively bound pairs IT231 l24l l25l l26ll . 
This picture is also supported by the pair probability pA,2(t) 
depicted in Fig. [6] the probability of finding two bosons in 
the same well remains close to one with small modulations, 
that can be interpreted as attempted single-particle tunneling. 
Although our effective Bose-Hubbard model (eq. (|6]l) is not 
strictly valid in this case anymore, it gives the right tendency 
for the tunneling times: T ~ 2-k8u/AJ 2 rs 260, compared to 
the numerically exact value T m 220. 
The dynamics of the heavier B bosons 

Also the dynamics of the B boson becomes more intriguing 
for an increasing mass ratio, as they begin to move due to the 
backaction of the lighter tunneling atoms. They remain fixed 
at the center of the trap in the limit of vanishing mass ratios 
a = 0.001 (Fig. [6j. For slightly increased a = 0.01, they 
carry out tiny oscillations around zero and only begin to move 
noticeably for larger mass ratios a sa 0.1. Their movement 
is always in reverse to the A atoms, which is visible com- 
paring their density evolutions. For instance, at a = 0.12, the 
lighter bosons tunnel as attractively bound pairs on a long time 
scale, hence also the dynamics of the B bosons takes place at 
long scales. This behavior is obvious considering the fact that 
they are strongly repelled by the A bosons and move in or- 
der to avoid them. By extension, this dynamics in reverse to 
the tunneling atoms is present for all intra-species interactions 
strengths gA- 
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Figure 8: (color online) Inter-species tunneling dynamics of Na = 2 
bosons with intra-species interaction g\ ~ 0.5 through a strongly re- 
pulsive B atom (t/AB = 8) and increasing mass ratios a — 0.001 
(— ), a = 0.01 (- - -), and a = 0.12 (••■); Top: Relative A popula- 
tion pA,ti(t) and Bottom: Pair probability pA,2(t) 



B. Intermediate interactions gA ~ 0.5 

The case of two moderately interacting bosons with gA = 
0.5 reveals an entirely different behavior as illustrated in 
Fig. [8] Near the static limit (a = 0.001), the atoms tunnel 
on very long time scales and only feature a minute faster os- 
cillation on top. This pattern is predicted by the spectrum in 
Fig. [7] two very different time scales determine the dynamical 
evolution. While the upper two lines Ei^ lie very close, lead- 
ing to a long tunneling time T = 2tt/lui2, the larger splitting 
of the lower lines woi only effects a tiny fast oscillation on top 
of it. At release time t = 0, both bosons are prepared in the 
same well and, due to their additional interaction energy, are 
off resonance with the configuration, where each well is oc- 
cupied by one boson. As the system is non-dissipative, no en- 
ergy can be lost, and it is never possible to access a state with 
only one boson in each well, although being energetically fa- 
vorable. Tunneling of single atoms is therefore strongly sup- 
pressed and only the tunneling of pairs is possible. The pah- 
correlation pA,2(t) confirms this picture: it roughly remains 
around 100 %, apart from small reductions due to attempted 
single-particle tunneling. This indicates that it is always very 
likely to find both A bosons together either on the left- or the 
right-hand side of the trap and underlines the stability of the 
pair. This behavior is reminiscent of the tunneling of repul- 
sively bound atoms pairs, that has been studied for external 
barriers, iflll^ 

If we slightly increase the mass ratio a = 0.01, the time 
evolution remains practically the same, as such small changes 
on the energy levels do not have a visible effect due to the very 
different tunnel frequencies wrji , (^12- 
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Figure 9: (color online) Inter-species tunneling dynamics of JVa — 2 
fermionized bosons (g\ = 25) through a strongly repulsive B atom 
(ffAB = 8) and increasing mass ratios a = 0.001 ( — ), ot = 0.01 (- - 
-), and a — 0.12 (• • ■ ); Left: Normalized A population pA,n(t) and 
Right: pair probability pa, 2(f) 



For much higher mass ratios a = 0.12, the repulsively 
bound pair finally breaks up and the sinusoidal oscillations 
of the A bosons are recovered. This behavior is predicted by 
the spectrum in Fig. [7] where w i ~ ^12 an d consequently 
the dynamics is governed by only one frequency. The effect 
becomes even more lucid from the angle of the pair proba- 
bility: instead of remaining around pa, 2 (t) ~ 1, it oscillates 
between one and one half. These properties are qualitatively 
predicted within the effective Bose-Hubbard model ©: the 
induced on-site attraction Su < « —0.34 roughly cancels 
the true intra-species repulsion uq > « 0.25 which leaves 
us with an effectively noninteracting system that carries out 
Rabi oscillations with a period T = tt / ps 65. 



C. The fermionization limit 

The dynamics becomes more involved as we increase the 
intra-species interaction up to the fermionization limit qa = 
25. Needless to say that in this strongly correlated limit the 
lowest band approximation as stated in eq. (0 breaks down 
and we have to find another way to explain the dynamics. 
According to the Bose-Fermi map jfj, hard-core bosons and 
fermions have an identical energy spectrum and we seek to un- 
derstand the tunneling dynamics of the bosons by effectively 
treating them as fermions. As a consequence of the Pauli prin- 
ciple, (spinless) fermions in one well cannot occupy the same 
band, and consequently populate the Aa lowest bands in the 
initial state. The probability of finding a particle at the right- 



hand side can therefore be written as 

, N A -1 

PA,r(*) = F I sin 2 (j(«i). (8) 

A 0=0 

As visualized in Fig.|9]for Aa = 2, the only two frequencies 
that contribute for a ~ are J*- ' and JW. Also the pair 
probability p\,2(t) supports the analogy between hard-core 
bosons and fermions: it takes nearly all values between one 
and zero, pointing out that the dynamics can be understood 
as the tunneling of two independent fermions with different 
frequencies. 

Recalling the effective Fermi-Hubbard model (eq. (0), 
with increasing a we expect an induced attraction between 
the tunneling fermionized bosons. In the intermediate regime 
(q = 0.01), this has only a small effect leading to a renor- 
malization of the contributing Rabi frequencies J 1 -^ with 
/3 = 0,1 for A A = 2. 

However, for a strongly increased mass ratios a = 0.12, re- 
markable effects emerge: single-particle tunneling is turned 
off-resonant and the two atoms move only on a time scale 
about four times longer than the Rabi oscillations. The tunnel 
dynamics closely resembles the one of a bound pair, a point 
of view that is supported by the pair probability PA,2(t) that 
stays remarkably close to unity, in marked contrast with the 
independent fermion tunneling for a w 0. In this light, it 
is tempting to think of the phenomenon as the tunneling of 
an attractively bound pair of identical fermions. The effec- 
tive Fermi-Hubbard model Q underpins this interpretation: 
the induced inter-band interaction <5u( 0101 ) is always negative, 
signifying that the fermionized bosons attract each other and 
only pair tunneling is resonant. While the effective model pre- 
dicts a tunneling period of T = 2n5v(° 10 ^ /U^ jU) « 320, 
the numerically exact result is T m 290. 

D. Higher atom numbers 

Although having concentrated on the special case of 
Aa = 2 light bosons so far, it would be interesting to 
consider higher atom numbers for several reasons. On the 
one hand, in a setup consis ting of a whole array of ID 
traps like described in J3, [E l32ll . number fluctuations may 
automatically admix states with Aa > 2. On the other 
hand we want to verify the validity of the effective Hubbard 
models for higher atom numbers. Furthermore, it would be 
interesting to connect the dynamics of our few-boson system 
to that of large atomic clouds. This shall first be discussed for 
an increased number of B bosons and thereafter for more A 
atoms. 

Higher numbers of B bosons: 
Increasing the number of B bosons does not change the tun- 
neling dynamics conceptually. In general, the tunneling times 
are enhanced and the B bosons are more inert in their reaction 
to the tunneling A bosons due to their increased total mass. In 
the static limit, a = 0.001, the heavier atoms display an effec- 
tive barrier for the lighter species, which is proportional to the 
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Figure 10: (color online) Inter-species tunneling dynamics of Na = 
2 bosons through 7V~b = 2 atoms for (jab = 8), g a =0 and increas- 
ing mass ratios: Normalized A population pA,n(t) for a = 0.001 
(— ), a = 0.02 (- - -), and a = 0.12 (• • • ) 

number of B bosons, rendering the tunnel coupling J between 
different sites weaker. The tunnel period T = it/ is ac- 
cordingly enhanced (see Fig. ITOt. For an increased mass ratio 
q = 0.01 (not shown), the beatings known from Ab = 1 are 
not observable: the total mass of the B bosons is larger, which 
results in a weaker backaction to the tunneling component. 
Hence the beatings can only be recovered for larger a = 0.02 
as displayed in Fig. \TU\ In the limit of larger mass ratios 
a = 0.12, we again observe tunneling of attractively bound 
pairs of A bosons on a longer time scale T ~ 27t<5m/4 j(°) 2 
due to the smaller tunnel coupling 
Higher numbers of A bosons: 

The weakly -interacting behavior of an increased number of 
A bosons does not differ conceptually from the case of lower 
particle numbers. In the limit of a — > 0, where the B bosons 
can be treated as static barrier for the A bosons, eq. (O carries 
over to 

1 ( NA \ 

PA,r(*) = ^ I 1 + X] c "i c "(*™l Q ( a; )l*n) COs(w„ TO t) J , 

where Aa modes have to be taken into account. For vanishing 
intra-species interactions gA = and a vanishing mass ratio 
a « 0, the same simple Rabi oscillations as for smaller atom 
numbers arise (see Fig.fTTTi. As we slightly increase the mass 
ratio to a = 0.01, the equidistance between the contributing 
levels is again lifted and beats evolve on top of the sinusoidal 
oscillations. 

As we move to larger mass ratios a = 0.12, the 
lowest-lying of the JVa states have merged into doublets, 
where the highest eigenstate pair is roughly described by 
|*a) = -^(\N L ,0 R ) ± \0 R ,N L )). This effect results 
from the induced intra-species attraction Su and can be 
extrapolated from the behavior for Aa = 2 in Fig. [7] Their 
energy splitting in this two-mode model has been determined 

to be lu oc pffrjr ( M ^-) NA HI- The tunnel P eriod 
is thus expected to grow exponentially with the number 
of particles as already foreshadowed in Fig. [TT] for the 
case of Na = 3 bosons. Extrapolating this behavior to 
higher particle numbers, we finally reach the condensate 
regime (Na 3> 1), where complete population transfer 




20 40 60 80 100 120 140 
time t 



DC 0.7 




20 40 60 80 100 120 140 
time t 




20 40 60 80 100 120 140 
time t 



Figure 1 1 : (color online) Inter-species tunneling dynamics of Na = 
3 bosons through a strongly repulsive B atom (jab = 8) and in- 
creasing mass ratios: Relative A population pA,rt(t) for a = 0.001 
(— ), a = 0.01 (- - -), and a = 0.12 (■ • • ), (a): g A = 0, Inset: 
Long-time behavior a = 0.12, (b): gx = 0.5, (c): gx = 25 



becomes completely inaccessible. The same prediction is 
made within the effective Bose-Hubbard model that yields a 
tunneling period of T w 2ttu 2 /3 J {0)3 ~ 2500 for three A 
bosons for a mass ratio a = 0.12, compared to the numeri- 
cally exact value T m 1800, depicted in the inset in Fig. fTTT a). 

The same scaling of the tunneling period can be observed 
for an intrinsically interacting system in the case of interme- 
diate coupling strengths gA ~ 0.5: for A^a = 3, we observe 
long tunneling times that increase exponentially with the 
number of A bosons. Only tiny corrections are remarkable 
for small increments of the mass ratio (a f=a 0.01). It is only 
for larger mass ratios a sw 0.12 that the induced attraction 
Su is sufficiently strong to balance the true repulsion uq, 
and we recover Rabi oscillations, already known from the 
case of Aa = 2 bosons in Fig. [TTT b). This behavior can be 
extrapolated to larger particle numbers, where the induced 
attraction 8u always cancels u at approximately a sa 0.12. 

For strong intra-species repulsions gA = 25, the fermion- 
ized bosons can again be mapped onto noninteracting 
fermions and eq. $Q yields the proper description. For Aa = 
3 hard-core bosons, three Rabi frequencies contribute, ex- 
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plaining the picture that emerges in the limit of vanishing 
mass ratios a w (see Fig. [TTT c)). In the intermediate 
regime (a = 0.01), beatings on each of the contributing 
"Rabi frequencies" evolve, which do not qualitatively 
alter the overall tunneling dynamics. For higher mass ratios 
a = 0.12, we can again observe the tunneling of attrac- 
tively bound pairs of fermions, which is properly accounted 
for in the effective Fermi-Hubbard model. Due to energy con- 
servation, only terms that conserve the number of atoms in 
each band contribute to the inter-species tunneling dynam- 
ics, such that the on-site interaction Hamiltonian reduces to 
H on - site = E s , Q</ 3'^ (Q/3Q/3) ">f- All contributions for 
N A = 3, namely (5v (0101) ,(5u (0202) ,(5w (1212) , are negative and 
lead to a tunneling time of T « 1600, compared to the numer- 
ically exact result T « 1150. 



V. CONCLUSION 

In summary, we have studied two-component mixtures of 
one-dimensional few-boson systems under harmonic confine- 
ment. In our setup, the two species experience very different 
confinement lengths, realized via different masses (in anal- 
ogy, this can also be achieved by a variation of the trapping 
frequencies). By varying the mass ratio, different regimes, 
from the complete localization of one species at the center 
of the trap to equal masses for both constituents, have been 
investigated. The inter-species interaction has been fixed at 
a large repulsive value, while the full crossover from weak 
intra-species coupling to the fermionization limit has been in- 
vestigated for the lighter bosons. 

In the limit where the two species experience a very different 
degree of localization, the heavier one can be considered as an 
effective barrier for the lighter species and we can map this to 
a system of identical bosons in a double well. On easing the 
localization requirement of the heavier bosons, an additional 
correlation between the lighter species arises. For this case, 
we rigorously derive an effective lattice model and show that 
one can account for the mobility of the barrier in terms of an 
induced attraction. For the ground state properties, this may 
lead to strong intra-species correlations even in the case of 
intrinsically noninteracting bosons, resulting in an increased 
probability of finding several atoms at the same position. On 
the other side, it is possible to cancel an initial intra-species 
repulsion by means of the induced attraction, leading us back 
to an effectively uncorrelated system. Also the tunneling dy- 
namics of the mobile component through the effective barrier 
reveals a similar physical behavior. In the limit of complete 
localization of the heavier species, the dynamics through the 
effective barrier bears a strong resemblance to tunneling in a 
double well. Increasing the mass ratio, it is possible to go from 
Rabi oscillations in a system with vanishing intra-species in- 
teractions, to a regime, where only tunneling of attractively 
bound pairs is possible. Notably, this is even realizable for 
fermionized bosons, that, despite of being strongly repulsive, 
begin to tunnel in attractively bound pairs through the heav- 
ier species. This dynamical behavior is not limited to systems 
with only one impurity and by extension, it would also be pos- 



sible to create and investigate an effective lattice of localized 
atoms, allowing for the study of highly-correlated quantum 
systems. 



Appendix A: COMPUTATIONAL METHOD 

The results described in this paper are based on the nu- 
merically exact Multi-Configuration Time-Dependent Hartree 
method ll34l [35||_that has been successfully applied to few 
boson systems 111 S US [U US [H HI [H. Its prin- 
cipal idea is to solve the many-body Schrtidinger equation 
i*&(t) = H^(t) as an initial-value problem by expanding the 
solution of the wave function in terms of direct (Hartree) prod- 
uct states $j(f ) = (t) ® .. ® <fj N (t) 



(Al) 



where the multiindex J = (ji,...,jjv) runs over all admis- 
sible configurations C = { J = (jx, Jjv)|1 < ji < s } with 
s denoting the maximum number of required basis functions. 
The single-particle functions {yjj(f)|l < j < s} that are used 
to built up the product states are still unknown and represented 
in a fixed primitive basis on a grid, which is in our case pro- 
vided by harmonic oscillator functions. 

Applying the Dirac-Frenkel variational principle yields 
equations of motion for both the coefficients Aj(t) and the 
single-particle functions (fj(t). Integrating this system of dif- 
ferential equations leads to the time evolution of the system 
via (IAU . In the above expansion not only the coefficients 
Aj(t), but also the single-particle functions (fj(t) are time 
dependent, which results in a basis set {$,/(£)} that is vari- 
ationally optimal at each time t and can thus be kept rela- 
tively small, lending this method its high efficiency. Still, 
the exponential growth of the numerical effort limits the ap- 
proach to particle numbers N < 10 depending on the number 
of single-particle functions necessary to describe the inter- 
particle correlations. The peculiarity of the system at hand, 
a two-component mixture of bosons, is the indistinguishabil- 
ity within each species, which manifests in identical single- 
particle functions within each subset Ka = {1, Na} and 
Kb = {Na + 1, N}. The bosonic permutation symmetry 
is not a priori included in the ansatz (IA 1 1 and is ensured by 
the correct symmetrization of the expansion coefficients Aj. 

Besides time-dependent dynamical investigations, this pa- 
per also contains stationary-state calculations, which can be 
obtained employing the so-called relaxation method II40I1 . The 
principal idea behind this approach is to propagate an ini- 
tial wave function ^'(O) by the non-unitary, imaginary time- 
propagator U(t) = e~ Hr . For r — > oo, this damps out 
any contribution from excited states by the exponential fac- 
tor e -( E M- E «) T and we are left with the ground state. This 
basic principle is very sensitive to numerical instabilities, and 
the scheme that one relies on in practice is the more robust 
improved relaxation [41]. The starting point of this ansatz 
is to minimize (\E'|7J|^I') with respect to both the coefficients 
A j and the orbitals ipj. The resulting eigenvalue problem 
is solved iteratively by first diagonalizing the Hamiltonian in 
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some fixed basis to obtain the coefficients Aj. The orbitals ipj 
are then optimized by propagation over a short period in imag- 
inary time and the Hamiltonian matrix is subsequently rebuilt 
with these new orbitals. This procedure is then repeated until 
the desired accuracy is obtained. 
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